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ABSTRACT 

After 3 years of quiescence, the globular cluster NGC 6440 exhibited a bright 
transient X-ray source turning on in August 2001, as noted with the RXTE All- 
Sky Monitor. We carried out a short target of opportunity observation with 
the Chandra X-ray Observatory and are able to associate the transient with the 
brightest of 24 X-ray sources detected during quiescence in July 2000 with Chan- 
dra. Furthermore, we securely identify the optical counterpart and determine 
that the 1998 X-ray outburst in NGC 6440 was from the same object. This is 
the first time that an optical counterpart to a transient in a globular cluster is 
securely identified. Since the transient is a type I X-ray burster, it is established 
that the compact accretor is a neutron star. Thus, this transient provides an ideal 
case to study the quiescent emission in the optical and X-ray of a transiently ac- 
creting neutron star while knowing the distance and reddening accurately. One 
model that fits the quiescent spectrum is an absorbed power law plus neutron 
star hydrogen atmosphere model. We find an intrinsic neutron star radius of 
nlill km and an unabsorbed bolometric luminosity for the neutron star atmo- 
sphere of (2.1 ± 0.8) X 10^^ erg s^^ which is consistent with predictions for a 
cooling neutron star. 

Subject headings: globular clusters: individual (NGC 6440) - stars: neutron - 
X-rays: binaries 

^Astronomical Institute, Utrecht University, P.O. Box 80 000, 3508 TA Utrecht, the Netherlands 
■^SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, the Netherlands 
■^email jeanz@sron.nl 

''^Center for Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139-4307 
^ Chandra fellow 



- 2 - 



1. Introduction 

Currently, thirteen bright X-ray sources are known to inhabit twelve Galactic globular 
clusters (for a recent overview, see Sidoli et al. 2001; the 13th source was recently discovered 
by White & Angelini 2001). They are low-mass X-ray binaries (LMXBs) of which at least 
11 harbor a neutron star. Four are transient in nature: those in Liller 1, Terzan 5 and 
6, and NGC 6440 (one could argue that Terzan 1 contains a transient source as well with 
an activity lasting decades, see Guainazzi, Parmar & Oosterbroek 1999). Thus far, none 
of these transients had a secTirely identified optical counterpart. The prospects for finding 
one were best for the transient in NGC 6440 because it has a relatively low reddening of 
E{B — V) = 1.07 and a non-extreme distance of 8.5 ± 0.4 kpc (Ortolani, Barbuy & Bica 
1994). The transient in NGC 6440 was seen in outburst in Dec. 1971 (Markert et al. 1975; 
Forman, Jones & Tananbaum 1976) and Aug. 1998 (In 't Zand et al. 1999; Verbunt et al. 
2000). From ROSAT HRI locahzations of two sources in quiescence and optical observations 
during the 1998 outburst, two candidate optical counterparts were identified, VI and V2 
(Verbunt et al. 2000). V2 was evaluated to be the most hkely one with B Ri 22.7 during 
outburst. 

Pooley et al. (2001) observed NGC 6440 for 23 ksec with the Chandra X-ray Observatory 
on July 4, 2000, when the transient was quiescent. They found 24 point sources with 0.5-2.5 
keV luminosities between 1.8x10^'^ and 9.9 x 10^^ erg s~^, which they designate CXI through 
CX24. The brightest source, CXI, coincides with the optical counterpart candidate V2. 

The transient in NGC 6440 erupted again in August 2001, as detected with the All-Sky 
Monitor (ASM) on the Rossi X-ray Timing Explorer RXTE^. We proceeded to secure the 
identification of the optical counterpart through swift, high-resolution X-ray observations, 
and report the results in the present paper. A preliminary report was already given by In 't 
Zand et al. (2001). 

2. Observations and source identification 

The ASM (Levine et al. 1996; Bradt et al. 2000) monitors NGC 6440 for 10 months 
of each year at a detection threshold of about 15 mCrab in the 1.5 to 12 keV band per day 

of observation. The transient exceeded this threshold during two weeks starting on August 
2, 2001. The peak flux was 60 ± 10 mCrab. Continued ASM monitoring showed a reflaring 
episode starting on Sep. 10, with fluxes up to 0.2 Crab and lasting for already almost 2 



^as announced by the ASM team at http://xte.init.edu 
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months at the time of this writing. The Wide Field Cameras (WFCs) on BeppoSAX (Jager 
et al. 1997; In 't Zand 2001) also observed NGC 6440 during this episode and detected at 
least 10 type-I X-ray bursts. As soon as the community was alerted of the outburst (on 
August 13), we proposed X-ray observations with Chandra and optical observations with the 
European Southern Observatory's New Technology Telescope in Director's Discretionary 
Time. 

The Chandra (Weiskopf, O'Dell & van Speybroek 1996) target-of-opportunity obser- 
vation (TOO) was approved and carried out on August 18.9 UT for an exposure time of 
2.6 kscc with the Advanced CCD Imaging Spectrometer (ACIS; Carmirc, Nousck & Bautz, 
in prep.) and no grating. Despite the anticipated heavy pile up, we chose this configuration, 
and not HRC, to optimize the sensitivity such that not only the transient but also a few 
of the other sources in NGC 6440 would be detected (ACIS has 2 to 4 times more effective 
area for 0.3-8 keV photons than HRC-I). Thus, the association of the transient with one of 
the quiescent sources would become independent from the Chandra attitude accuracy. Only 
a 2' X 8' subarray of the back-illuminated S3 chip was read so that the CCD frame time 
was 0.8 s instead of 3.2 s and pile up effects were reduced. Figure 1 shows the obtained 
image, zoomed in on the globular cluster and compared to the observation during quiescence 
(Pooley et al. 2001). It is evident that the transient was active. Furthermore, it is closest to 
CXI, the brightest of the sources detected during quiescence. Based on the Chandra aspect 
solution for the 2001 observation (accurate to about 0.5") and the astrometry carried out 
by Pooley et al. (2001) for the 2000 observation, the angular distance between the transient 
and CXI is 0.14" (the position of the transient is determined through weighing pixel position 
with pixel values). Of the remaining sources, CX2 and CX5 are clearly detected during the 
TOO as well, with photon rates in the full bandpass of (8 ± 2) x 10~^ and (6 ± 1) x 10~^ s~^. 
These rates compare well with those during the earlier observation ((7.6 ± 0.6) x 10~^ and 
(4.3 ±0.4) X 10^^ s^^ respectively). The angular distances of the transient to CX2 and CX5 
are consistent with that to CXI in the July 2000 image within 0.28±0.13". The distance is 
inconsistent with all of the other 23 sources from the July 2000 observation by at least 15a. 
This confirms the transient-CXl association. 

A total of 2.6 X 10^ events was detected within 8.0" from the source centroid. In 
order to estimate the flux and avoid piled-up pixels, we extracted a spectrum from the 
wings of the point spread function (PSF) for pixels that are between 2.5 and 8.0" from 
the centroid. We ignore negligible contributions from other sources in this annulus. An 
absorbed power law fits the 0.7-8.0 keV spectral data well, with a photon index of 1.2 ± 0.1, 
A^H = (8.2 ±1.6) X 10^^ cm~^ and xt = 1-40 (i^ = 75 degrees of freedom). After a correction 
for the fraction of the PSF covered (9%), the flux is 1.1 x 10~^° erg s~^cm^^ (0.5-10 keV) and 
the luminosity 9 x 10^^ erg s~^ (1.1 x 10^^ erg s~^ unabsorbed). The statistical uncertainty 
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is 7%, that due to PSF uncertainties an estimated 10%. The flux is roughly 10 times fainter 
than at the peak and 10^ times brighter than during quiescence (see below). The source 
shows no strong variability, with an upper limit of 20% at 90% confidence for 100 s time 
scales. No narrow features are apparent in the power density spectrum. 

The position of CXI (Pooley et al. 2001) is consistent with that of the candidate optical 
counterpart V2 in Verbunt et al. (2000). They are 0.46" apart while the uncertainty in the 
optical position alone is 0.5". We conclude that V2 is the secure optical counterpart to the 
transient in NGC 6440, that CXI is the quiescent X-ray counterpart, and that the 1998 and 
2001 outbursts were from the same object and not from any of the other 3 or 4 transient 
low-mass X-ray binaries in quiescence that are likely to be present in the cluster (Pooley 
et al. 2001). Most importantly, since type-I X-ray bursts were detected from the transient 
during both 1998 and 2001 outbursts, it is now flrmly established that CXI is a quiescent 
neutron star. 

NGC 6440 was observed for us in Service Time with the New Technology Telescope 
(NTT) at La Silla 3.5 d after the Chandra observation using the Superb Seeing Imager 
SUSI2. The seeing was around 1". Images were obtained in the Bessell R and B fllters. The 
analysis of these data is complicated by proper motions of a number of stars in the fleld. 
However, the raw images do show V2 to be active again in B at a level comparable to the 
1998 outburst. 



3. The quiescent X-ray spectrum 

Quiescent neutron stars (qNSs) exhibit X-ray emission at typical luminosities of 10^^ — 
g-i (g g^ Asai et al. 1998, McClintock 2001). A plausible explanation for this 
emission is thermal radiation from the atmosphere of a neutron star that is cooling from 
the previous outburst (Brown, Bildsten & Rutledge 1998). We extracted the spectrum of 
the transient during quiescence by accumulating all 251 photons within 1.5" from the source 
centroid and applying background corrections as found elsewhere in the image (amounting 
merely to a 2.3 photons subtraction). To ensure applicability of the statistic, we rebinned 
the data so that each channel contains at least 20 photons except for the highest-energy 
bin which contains 15 photons. We used CIAO 2.1.3 Science Threads for the background 
subtraction and response matrix generation. The spectrum was fltted between 0.3 and 6.7 
keV (only 2 photons were detected above this range) with a model combining a power law 
and a neutron star hydrogen atmosphere model as calculated for the non-magnetic case 
by Zavlin, Pavlov & Shibanov (1996). Standard absorption by cold interstellar matter of 
cosmological abundances was applied following Morrison & McCammon (1983). Due to 
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the limited statistical quality, the fit errors for the atmosphere temperature and radius are 
reasonably constrained only if the mass of the neutron star, the distance, the interstellar 
hydrogen column density A^h and the power law photon index are fixed. We did so at values of 
1.4 Mq (for a 'canonical' neutron star), 8.5 kpc (this is well determined), A^h = 8.2x10^^ cm~^ 
and 1.44 (these two values were determined from the 1998 outburst; In 't Zand et al. 1999), 
respectively. A^h has an uncertainty of roughly 40%. The only three parameters left free 
are the neutron star radius, the atmosphere temperature and the power law normalization. 
The resulting fit is good, with xl — 0.758 (for u — 9 degrees of freedom). The neutron 
star parameters are strongly dependent. The fitted temperature is 92^^33^ eV and the radius 
ITtfa km (90% confidence). The radius to an observer at infinity is 19l|^ km. A graph 
of the spectrum is presented in figure 2 while in figure 3 the 90% confidence region in 
temperature/radius space is presented. If A^h is left free, it converges to 8.2^1'^ x 10^^ cm~^ 
and the range of allowed NS radii surpasses the 100 km mark. This upper radius margin is 
also quite sensitive to the power law index which is poorly constrained by the data. 

The inferred radius is consistent with that expected for a neutron star (e.g., Lattimer 
& Prakash 2001). We note that the model implies an intrinsic radius that is 20% larger 
if the NS mass is 1.8 instead of 1.4 Mq. Fits with some other models are nearly as good 
(table 1), so the data are not discriminate in that sense. The absorbed 0.5-7 keV flux is 
(5.9 ± 0.4) X lO-^'' erg s-^cm-^, the unabsorbed flux (2.0 ± 0.6) x 10"^^ erg s-^cm-^ (the 
large error reflects the uncertainty in A^h)- The power law contribution in the latter flux 
is 20±6% (90% confldence); this is only 8% in 0.5 to 2 keV. The bolometric unabsorbed 
luminosity of the NS atmosphere is (2.1 ± 0.8) x 10^^ erg s"^ (the contribution from the 
uncertainty in the distance is small). 



4. Discussion 

Pooley et al. (2001) flnd that there are 4 or 5 likely quiescent LMXB transients in 
NGC 6440. Similar conclusions are currently being reached from Chandra observations of 
other globular clusters such as 47 Tuc (Grindlay et al. 2001a), Liller 1 (Homer et al. 2001), 

NGC 6652 (Hcinkc, Edmonds & Grindlay 2001), u Ccn (Rutlcdgc et al. 2001) and NGC 6397 
(Grindlay et al. 2001b). Furthermore, for the first time two bright low-mass X-ray binaries 
have been identified in a single cluster (M15; White & Angelini 2001). Our observations 
show unequivocally that the second and third transient outburst in NGC 6440 (in 1998 and 
2001) are from the same object. Due to a limited positional accuracy of the source of the 
1971 outburst, it is not clear whether the same source was responsible for the 1971 and 
the 1998/2001 outbursts. At first, indirect evidence suggested that the source is different. 
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because the duration of 1.5 to 3 months and the peak flux of 250 mCrab deviate substantially 
from the values for the 1998 outburst and for the initial phase of the 2001 outburst. However, 
at the time of this writing the 2001 outburst shows activity for already 3 months and a peak 
flux in excess of about 0.2 Crab units, according to ASM and WFC measurements. 

The recurrence time between the two most recent outbursts is unexpectedly short. This 
may mean that the source has become more active, or alternatively that earlier outbursts 
have been missed. It is not uncommon for a previously quiet source, to become active in a 
series of transient outbursts (e.g. GRO J1744-28 and GRO J1655-40). However, the 1998 
outburst of the transient in NGC6440 was so faint and short, that similar outbursts may 
easily have passed undetected. We note that the RXTE ASM and BeppoSAX WFCs became 
operational only in 1996. Similar coverage as the ASM in the future should resolve this issue. 

This is the first time that an optical counterpart has firmly been established for an X-ray 
transient in a globular cluster. Furthermore, it is the first time that a quiescent neutron star 
has been identified without source confusion in a globular cluster. These identifications are 
important because they present the case of a quiescent neutron star transient in the optical 
and X-rays for which an accurate distance is known. 

The spectral study of quiescent neutron stars is foremost interesting, as Brown et al. 
(1998) pointed out, because it may result in an accurate enough neutron star radius deter- 
mination to be able to constrain the equation of state of matter in the star's core. Measuring 
the radius to an accuracy of about 0.5 km can exclude 50% of the proposed equations of 
states (Lattimer & Prakash 2001). For neutron stars in globular clusters as compared to 
those in the Galactic disk, one source of uncertainty is almost eliminated — the distance. 
The distance to NGC 6440 is accurate to 5%. If other uncertainties (in the spectral pa- 
rameters) had been at least an order of magnitude lower, the prospects for constraining the 
equation of state for the transient in NGC 6440 would have been good. However, this is not 
the case. This qNS would have to be observed for a few million seconds with Chandra before 
it becomes worthwhile to consider this analysis, and that is without considering uncertainties 
in the spectral models. 

Given a peak bolometric luminosity between 4 x 10^ and 1.6 x 10^^ erg s~^ (or 50 to 
200 mCrab in 2-10 keV given the broad-band spectrum measured by In 't Zand et al. 1999), 
an e-folding decay time between 5 d (Verbunt et al. 2000) and 30 d, and assuming that 
during accretion 100% of the gravitational potential energy is released through radiation, 
the recurrence time of 2.95 yr imphes a time-averaged accretion rate of < M >Ki 2 x 10~^^ 
to 4 X 10^^° Mq yr~^. Following Brown et al. (1998), the quiescent time-averaged luminosity 
due to thermal emission from the cooling NS is predicted to be between 1 x 10'^'^ and 2 x 
-|^g34 gj,g g-1 Q^j, measurement of 2.1 x 10^^ erg s~^ is consistent with this rough prediction. 
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also if less than 100% of the hberated potential energy is released as radiation. 
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Table 1: Goodness of fit for a variety of spectral models to the 12-cliaiinel data between 0.3 
and 6.7 keV. 



Modelt 


free 








params. 






wa*nsa 


T,i? 


2.74 


10 


wa*(nsa+po) 


T, i?,po-norm 


0.758 


9 


wa*bb 


T,bb-norm 


3.00 


10 


wa*(bb+po) 


T,bb-norm,po-norm 


0.987 


9 


wa*diskbb 


A^H,7in,norm 


1.94 


9 


wa*diskbb 


Tin, norm 


2.44 


10 



Wa = interstellar absorption following Morrison & McCammon (1983), nsa = NS H- 
atmosphere emission according to Zavlin et al. 1996, po = power law, bb = black body, 
diskbb = disk black body according (Mitsuda et al. 1984) 



Fig. 1. — Chandra ACIS-S images of NGC 6440 during the transient active state (left; 2.6 
ks exposure) and quiescence (right; 23 ks exposure). The grids are fixed to the equatorial 
coordinate system and defined in the same manner for both cases. Pixels are O'.'S on the 
side. The left panel shows the bright transient with its large point spread function. The 
hole in the middle of the PSF and the streak are artifacts introduced by the brightness of 
the transient. The bright transient is coincident with quiescent source 'CXI' (Pooley et al. 
2001). CXI is the brightest X-ray source in the globular cluster during quiescence with a 
0.5-7 keV unabsorbed luminosity of (1.5 ± 0.5) x 10^^ erg s~^. 
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Fig. 2. — The spectrum of the quiescent emission of the transient, as obtained in July 2000. 
The crosses refer to the measurements, with the vertical lines indicating the errors, and the 
histogram refer to the fluxes per channel as predicted by the best-fit NS atmosphere model 
plus power law, after resetting the power law contribution to zero to illustrate its impact 
at high energies. The resetting also illustrates that both spectral components cover fairly 
distinct photon energies. 
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Fig. 3. — Contour plots of 90%-confidence (R,T) region of NS atmosphere model fit to CXI 

spectrum during quiescence between 0.3 and 6.7 keV, for an observer at infinity (left) or at 
the NS surface (right). It is assumed that the distance is 8.5 kpc and A^h = 8.2 x 10^^ cm~^. 
The location of the contour is independent of NS mass for an observer at infinity. 



